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Abstract 

We present a cosmological model containing a cosmological constant A and a 
component with an inhomogeneous equation of state. We study the form of the 
inhomogeneous equation of state for which the model exhibits the relaxation of 
the cosmological constant, i.e. it asymptotically tends to the de Sitter regime 
characterized by a small positive effective cosmological constant. The effect of the 
relaxation of the cosmological constant is observed both for negative and positive 
values of A and for a range of model parameters. A special emphasis is put on 
the study of the details of the CC relaxation mechanism and the robustness of 
the mechanism to the variation of model parameters. It is found that within the 
studied model the effective cosmological constant at large scale factor values is 
small because the absolute value of the real cosmological constant is large. 



1 Introduction 

The understanding of cosmology has undergone a phase of intensive development 
during the last decade. It has been largely propelled by the arrival of more precise 
and abundant observational data [TJ [3J. The picture of universe implied by 
the new data, however, revealed that the unknown part of the composition of the 
universe is far greater that previously believed. It is interesting to notice that 
the progress achieved in a great deal consists in establishing the fact that we 
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know much less about our universe than we thought before the advent of the new 
observational data. The term unknown composition of the universe should not be 
understood too literally. It might really be the case that the universe contains 
a new physical component (or several of them), but it is also possible that the 
perceived unknown composition is just a manifestation of additional dimesions or 
modifications of interactions such as gravity. 

One of the most striking consequences of the unknown part of the composition 
of the universe is its present accelerated expansion which started at redshift values 
of the order 1. The accelerated expansion of the present universe is strongly con- 
firmed by the observational data. Presently available data, however, provide much 
less information on the actual cause of the acceleration. As already stated, the 
acceleration might happen owing to the existence of an unknown component with 
the negative pressure, referred to as dark energy (DE), or might be a consequence 
of the fact that our universe has a number of macroscopic dimensions different 
than 4 or that the laws of gravitational interaction are modified at cosmological 
scales. These possibilities represent some of the most studied options leading to 
the accelerated expansion, but they certainly do not exhaust the list of proposals 
for the explanation of the accelerated expansion. The concept of dark energy seems 
especially useful in the modelling of the acceleration mechanisms. Namely, even if 
the acceleration is not due to some physical component with a negative pressure, 
the framework of dark energy can be used as a very good effective description of 
the alternative acceleration mechanisms which especially facilitates the compar- 
isons of different approaches to the explanation of the accelerated expansion of the 
universe. 

It is also important to stress that the unknown composition of the universe is not 
entirely connected to the accelerated expansion of the universe. The observational 
data imply that our universe contains a significant component of nonrelativistic 
matter, also called dark matter (DM) which is important for the explanation of 
the growth and formation of the structures such as galaxies and clusters of galaxies 
that we observe in the universe today. The nature of dark matter also has not yet 
been firmly established. 

The dark energy component is primarily characterized by its negative pressure. 
A very large number of DE models have been proposed lately, of various degree of 
complexity, predictive potential and connection to fundamental physical theories 
[3]. The observational data still provide a lot of space for dynamical DE models, 
but the central place of the allowed parametric space is occupied by a very simple 
DE model, a so called ACDM model. This model assumes that the DE component 
is actually a small positive cosmological constant (CC). The cosmological constant 
is a well known concept present in the theory of general relativity (GR) from the 
very first years of the development of the theory [I3E1E!- An important observation 
is that the GR allows the existence of the CC, but it does not determine its size. So, 
in any cosmological model based on GR we are not concerned with the question 
whether the CC should be there or not, but with the problem of its size and 
sign. Since, therefore, the CC should already be an ingredient of cosmological 
model based on GR, it is very convenient to use such an object as a source of 
the acceleration of the universe. Indeed, an assumption of an existence of a small 
positive CC, together with the existence of dark matter, fits the observational data 
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very well. From the observational side things work well: we have a model with few 
parameters that uses familiar concepts and fits the data well. 

Fundamental quantum physical theories, however, predict various contributions 
to the observed value of the cosmological constant. In quantum field theory (QFT) 
there exist very large zero-point energy contributions for each of the quantum 
fields. There are also contributions from condensates such as Higgs condensate or 
QCD condensates. These contributions should be added to, in principle arbitrary 
value of CC allowed in GR. Any attempt of calculation of contributions to the 
cosmological constant reveals that their size is many orders of magnitude larger 
than the observed value of the CC. The number of orders of magnitude differs 
with the choice of effective QFT cutoff, but in any case we have differences which 
raise a lot of suspicion, to put it mildly. But, in the end, it is not that individual 
contributions matter, but their sum. In principle, for the classical contribution to 
the CC we can choose the needed value and reproduce the observed value of A. 
The problem is that all contributions have to cancel to very many decimal places 
for this mechanism to be effective. This problem, referred to as fine-tuning plagues 
this explanation of the observed value of the cosmological constant. It is sometimes 
also called "the old CC problem" . 

If theoretical considerations reveal such difficulties for the cosmological constant 
as a DE candidate, maybe we should opt for some of dynamical DE models. Even 
should the future observational data prove that dark energy (as a true component or 
as an effective representation of some other acceleration mechanism) is dynamical, 
it only relegates the CC problem to another level. Then we have to understand 
why the size of the CC is much smaller that the observed DE energy density or, 
possibly, why it is zero. The possibility that the CC is exactly zero would open the 
way to the solution of the CC problem by invoking some new symmetry. However, 
presently there is no proof for the existence of such a symmetry. 

Therefore, the CC problem is difficult and it goes even beyond the issue of dark 
energy. Indeed, apart from the drastic problem of the size of the observed CC, there 
is another problem related to the cosmological constant. Namely, presently avail- 
able observational data imply that the energy densities of nonrelativistic matter 
and the cosmological constant are of the same order of magnitude at present epoch 
of the expansion of the universe. The energy density of nonrelativistic matter 
scales very differently with the expansion than the CC. Therefore it is quite pe- 
culiar that these energy densities that have been very different in size in the past 
and will presumably be very different in the future (in the context of the ACDM 
model), are presently of comparable size. This problem is also called "coincidence 
problem". In this paper we shall mainly deal with "the old CC problem" whereas 
the "coincidence problem" will be only remotely commented. 

This paper further elaborates the mechanism of the relaxation of the cosmo- 
logical constant proposed in [8]. The relaxation of the cosmological constant is 
defined as a dynamical solution of the "old CC problem" without the fine-tuning 
of the parameters of the model. Essentially, in [8] we model the dark energy sector 
and study the asymptotic behavior at large scale factor values. The relaxation of 
the cosmological constant corresponds to the asymptotic de Sitter regime with a 
small positive effective cosmological constant. The mechanism is implemented in 
the framework of a cosmic component with an inhomogeneous equation of state 
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[9]. In this paper we expand the model of [8], study some of its limitations and 
examine the robustness of the CC relaxation mechanism. 



2 The cosmological constant relaxation model 

We consider a two component cosmological model containing a cosmological con- 
stant with the energy density p\ and an additional cosmological component with 
the energy density p. It is assumed that the universe is spatially flat, k = 0. The 
expansion of the universe is defined by the Friedmann equation 

H 2 = —(p A + p). (1) 

The evolution of the second component with the expansion of the universe is 
defined by the standard equation of continuity 

dp = -3(p + p)-, (2) 
a 

where its equation of state (EOS) has a nonstandard form 

p = wp-3( (H 2 + pr. (3) 

Here we take Co > and a is an arbitrary real parameter. The form of EOS given 
in ([3]) fits into the framework of the inhomogeneous equation of state [9]. The 
concept of the inhomogeneous equation of state was also studied in [TOj [HI E21 CG] ■ 
The role in the inhomogeneous DE equation of state in the process of structure 
formation was examined in [14j . In the remainder of this paper we study how the 
inhomogeneous equation of state contributes to the relaxation of the cosmological 
constant. The reference [9] (see especially the Appendix) shows that a possible way 
to understand the inhomogeneous equation of state is as an effective description of 
the modified gravity theories or time-dependent nonlinear viscosity. 

The modified theories of gravity study the extension of GR as a possible source 
of the acceleration mechanism active at present cosmological era [151 [T6] . An ex- 
ample of f(R) gravity, free of instabilities, [15l[16j[T7] was presented in [8] showing 
that the mechanism of the CC relaxation could be realized directly in f(R) modified 
gravity theories. A systematic study of modified gravity theories consistent with 
the solar system precision gravity tests [IS] is needed to establish the robustness 
of the CC relaxation mechanisms within the modified gravity theories. 

The concept of bulk viscosity, as a dissipation mechanism of imperfect cosmic 
fluid consistent with the symmetries of the FRW universe, was used for the study of 
various phenomena in cosmology \19 \ \2Q \ \21\. A very interesting possibility is that 
the bulk viscosity could potentially account for the present acceleration of the uni- 
verse without the presence of dark energy [22\ [23] . Here we consider generalization 
of the phenomenon of bulk viscosity. Namely, H is not a variable of state of the 
imperfect fluid and a general dependence of the fluid pressure on H represents a 
step out of standard framework of bulk viscosity. A more appropriate name would 
be time-dependent or nonlinear viscosity. This generalized concept of viscosity, 
however, proves to be very useful in the study of the present accelerated cosmic 
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expansion [23], peculiar properties of dark energy, including the phenomenon of the 
CC boundary crossing [25] as well as other interesting phenomena [26^27] . Further- 
more, it is of particular interest to investigate how the concept of viscosity combines 
with the concepts of braneworlds and modified gravity [281 EHl EOJ EH [32j [33] . 

The expressions ([T]), © and © can be easily combined to obtain a dynamical 
equation for the evolution of the Hubble function 

dH 2 + 3(1 + u,)^ (V - - ^(H 2 + P)A = . (4) 

a \ 3 1 + w J 

Further we scale all quantities of interest and introduce the following notation 



h = (H/H x ) 2 , s = a/a x , X = 8irGp A /3H x , £ = %irG^H% a 1] /{l+w), b = (3/H 2 x . 

(5) 

Here H(a x ) = H x . it is important to state that the value of a x can in principle 
take any value. It is not intrinsically constrained within the present model. Using 
this notation we can rewrite (Jl|) as 

s^ + 3(l + w)(h-\-£(h + b) a ) = 0, (6) 
as 

with h(l) = 1 defining the initial condition. 

A thorough analysis of the model for the case 6 = was performed in [S] 0. 
There it was shown that the CC relaxation mechanism was active for a < 0. Here 
we proceed with the analysis of the same interval for a and examine the effects of 
the nonzero values of the parameter b. 

Next we consider the value a = —1 as a representative and an analytically 
tractable case. The equation (jHJ) now reads 

{k + h)dh - -3(1 (7) 



(h — h*x)(h — h*2) s 

Here h*i and h*2 stand for the zeros of the denominator of the expression at the 
left hand side of ([7]). Their respective values are given by the following expressions: 



1 

ti*i — 

and 



Ki = - (A-6 + V(A + 6) 2 + 4e) (8) 



^ 2 = ^(a-6- v /(A + 6) 2 + 4^ . (9) 

The differential equation (|7|) can be easily integrated and we arrive at the 
closed form solution for the dynamics of the scaled Huuble parameter with the 
scale factor: 

Here A\ = (b + ft*i)/(/t*i — K 2 ) and A 2 = —{b + /i*2)/(^*i — ^1*2)- 



1 Note that in [5] it was convenient to choose the parameter a somewhat differently that in the present 
paper. 
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Figure 1: The dependence of h on the scale factor a for the negative cosmological constant 
in the regime a < 0. The parameter values used are b = 0, a = — 1, A = —2000, £ = 0.02 
and w = —0.8. 



Before the analysis of the results presented above, we make a short summary of 
the results of paper [8] which correspond to a specific value 6 = 0. This is a starting 
point of our analysis since in this paper we are interested in how the nonvanishing 
value of parameter b modifies the CC relaxation mechanism observed in [8]. The 
dynamics of the Hubble function h depends on all model parameters a, A, £ and 
w. 

The Hubble function h as a function of the scale factor for a case of negative A 
with a large absolute value and other representative parameter values (£ > 0, a < 
and w > — 1) is given in Fig. [TJ For these intervals of parameters the dynamics of h 
is characterized by a very abrupt transition between a phase of expansion at small 
scale factor values where h ~ a~ 3 ( 1+w ) and a de Sitter phase at large scale factor 
values characterized by a small effective positive CC with h asym ~ ^-eff ~ £/|A|. 
With other parameters fixed, the scale of h before the transition grows with |A| 
and the value of ^- e ff decreases with |A|. The dynamics of h before the transition 
is not affected by the size of £ whereas its asymptotic value at large a grows with £. 
The choice of exponent a does not affect the behavior at small values of the scale 
factor, whereas the asymptotic value of h at large a decreases as a becomes more 
negative. Finally, as already stated, the value of w affects the behavior before the 
transition and the large a behavior of h does not depend on w. For other choices 
of parameters the model may exhibit other interesting types of dynamics which 
however do not correspond to the CC relaxation mechanism. 

The dynamics of the Hubble function h for a large positive A and representative 
parameter values (£ < 0, a < and w < —1) is given in Fig. [2l In this regime, 
both for small and large scale factor values we find de Sitter regimes, h ~ A at small 
a and h ~ — £/A at large a. These asymptotic regimes are again interconnected 
by an abrupt transition. The scale of the de Sitter regime preceding the transition 
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Figure 2: The dynamics of h as a function of the scale factor a for a positive A and 
a < 0. The values of the used parameters are b = 0, a = —1, A = 2000, £ = —0.02 and 



grows and the scale of de Sitter regime following the transition decreases with the 
size of A. The value of a does not affect the behavior at small a values and the 
scale of the de Sitter regime at large a grows as a becomes more negative. The 
asymptotic value of h grows with the size of |£| at large a whereas the dynamics at 
small a is not sensitive to the value of £. The value of parameter w does not affect 
the dynamics of h at large a whereas the approach to de Sitter regime at small a 
is sensitive to w. As in the case of negative A, for other parameter intervals the 
behavior of the model is different and the CC relaxation mechanism is not effective. 
It is also important to stress that in the case of positive A, the energy density p 
should be negative. This is a strong argument to consider the second component 
as an effective description of some other fundamental mechanism. 

Next we turn to the case of nonvanishing b. Since the values ([8]) and (J9|) deter- 
mine the asymptotic behavior of the model, let us further study their dependence 
on the parameter b. We generally assume that the parameter |A| is by far the 
largest parameter of the model. More precisely, we suppose that |A| 2 > |£| and 
|A| S> \b\. These assumptions allow us to make an expansion of the square root 
terms in (|8|) and ([9]) with the following results: 



The expressions differ for the cases of positive and negative A. For A > we 



w = 



-1.2. 







obtain 
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b 2 +4f 

h, 2 = -b- — . (14) 



4|A| 



On the other hand, for A < we have 



b 2 + 46 

/M = -6+^, (15) 
h*2 = A -rr-j— . (16) 



From the expressions (|13|) to (|16|) we can determine the effect of the parameter b 
on the asymptotic values. 

For a positive A, at very small values of \b\ when \b\ <C |£/A| we have h#i ~ A 
and /i*2 — — 6/A. On the other hand, for a sufficiently large \b\, where |6| 3> |6/A|, 
we have h*i ~ A and h*2 — —b- The asymptotic behavior at large scale factor 
values is determined by the value h*2 an d we see that at sufficiently large values of 
|6|, in the sense defined above, the parameter b determines the asymptotic behavior 
of the Hubble function h. It is important to notice that in this case de Sitter regime 
at large values of the scale factor is realized only for negative values of b. 

For negative values of A at small values for the parameter b, with |6| <C |£/A| 
we obtain and h*i ~ £/|A| and /t*2 — ^- For \b\ 3> |£/A| we further have h*i ~ —b 
and h*2 — A. These results again show that for a sufficiently large value of |6|, this 
parameter determines h*i which in turn controls the asymptotic dynamics of the 
Hubble function at large scale factor values. Again, to have a de Sitter regime at 
large scale factor values b has to be negative. 

A more careful analysis of the model immediately shows that for h — > —b the 
pressure of the component with the inhomogeneous EOS diverges. The dynamics 
of the model reveals that this singular point is never reached. The expressions © 
and Q show that the point where p would diverge is never reached during the 
evolution of the model. For a large |A|, the dynamics of h stabilizes at a value 
slightly above —b. 



3 Discussion 

The principal role of the considerations given above and the plots in Figures [3] and 
H] is to gauge the role of the parameter b in the model defined by dH)-®. Our 
primary goal is to find out if and in which extent does the finite value of b change 
the behavior of the model compared to the previously studied case corresponding 
to b = 0. The qualitative behavior of the model retains the pattern observed for 
the vanishing value of b: both for positive and negative values of the scaled CC 
parameter A there is a distinct and abrupt transition from the expansion at a high 
energy density to the de Sitter regime. Therefore, this specific signature of the 
CC relaxation mechanism is not lost with the addition of the additional parameter 
b. The asymptotic value of H 2 , however, depends on the interplay of all model 
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Figure 3: The behavior of h as a function of the scale factor a for different values of 
the parameter b. The used parameter values are a = —1, A = —1000, £ = 0.01 and 
w = —0.9. In this parameter regime, the value of b controls the asymptotic behavior of 
h at large a. 
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Figure 4: The dependence of the dynamics of h as a function of the scale factor on 
parameter b. The values of the parameters used are a = — 1, A = 1000, £ = —0.01 and 
w = —1.1. The asymptotic behavior of h at large a is determined by the size of b. 
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parameters. For a sufficiently small value of b, the asymptotic value is determined 
by the ratio of parameters £ and A (h = |£/A|). As b grows, the asymptotic value 
becomes fully dominated by the value of parameter b. 

Our aim is to look into a mechanism of the CC relaxation without fine-tuning. 
The case of vanishing b possesses certain appeal since there the effective positive 
CC is small because A is large in absolute value and also the parameter £ is not 
expected to be large. For parameter values which we would expect based on the 
fundamental theories the expected value of the effective CC at large scale factor 
values is a small positive number. There is no need for the parameters to be 
fine-tuned. 

For a large absolute value of A, sufficiently large b determines the asymptotic 
behavior of H 2 . If we eventually aim at explaining the observed value of the 
cosmological constant, the value of b should be small. Since there is no clear 
reason why the value of b should be so small, some form of fine-tuning reenters 
into the model. We have to introduce a small parameter b just to match the value 
of observed A e jy. Still, it is very important to stress that this value is very different 
from the real value A. Although a very large A is present in the model, it does 
not determine decisively the asymptotic behavior of the system. Furthermore, the 
parameter b plays the main role only because A is very large in absolute value. In 
a way the spirit of the CC relaxation mechanism is preserved: A universe with 
a large A finally tends to a de Sitter state characterized by a small A e ff. The 
principal difference to the b = case is that for a sufficiently large \b\ there is no 
strong argument why A e jj should be small. 

The model of this paper represents and extension of the model studied in [8], 
but it is still just a starting point towards a realistic cosmological model with the 
resolved CC problem. Essentially, both the model of this paper and [8] model the 
dark energy sector of the universe. Clearly, other components such as radiation 
and matter have to be added to create a realistic cosmological model and reproduce 
the standard eras of the evolution of the universe such as radiation dominated and 
matter dominated eras. The abruptness of the observed transition might pose a 
significant challenge to the construction of such a complete cosmological model. 
However, even in a model which contains the matter and radiation components, 
the presented CC relaxation mechanism should be efficient asymptotically since 
the energy densities of these components decay quickly with the expansion. Other 
important issues for further work are the timing of the transition, possible links 
to inflation and the growth of inhomogeneities in a universe in which the CC 
relaxation mechanism is active. 

As discussed in the Introduction, there is another important problem related 
to the size of the dark energy (or CC) density. The coincidence problem is not 
directly addressed in the present model since the matter and radiation components 
are not present in the model. Only in the model with all relevant components this 
issue could be addressed properly. 

Another very important issue is the physical motivation for the EOS given 
in ([3]). As already stated, the motivation could come from several directions of 
which we singled out modified gravity and generalized nonlinear viscosity. The 
elaboration of these topics might prove essential for a microscopic foundation of 
the mechanism exhibited by our model. 
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4 Conclusions 



The extension of the original model of the CC relaxation [8] presented in this paper 
allows us the study of the limits of the CC relaxation mechanism. This is achieved 
through the introduction of the new parameter b. For a sufficiently large \b\, we 
no longer have an explanation of the smallness of A e yj without fine-tuning. Still, 
even for larger values of |6| the asymptotic value is not determined by a large A, 
but some other small parameter, in particular b. In this sense, the spirit of the CC 
relaxation mechanism persists even for larger values of \b\. These conclusions show 
that the concept of the CC relaxation mechanism is robust, although for small 
or vanishing value of b the solution of the CC problem is more natural. These 
findings further support the study of other types of inhomogeneous EOS as a road 
to a complete cosmological model in which the CC problem is naturally solved. 

Acknowledgements. The author would like to thank B. Guberina for use- 
ful comments on the manuscript. This work was supported by the Ministry of 
Education, Science and Sports of the Republic of Croatia under the contract No. 
098-0982930-2864. 



References 

[1] A.G. Riess et al., Astron. J. 116 (1998) 1009; S. Perlmutter et al., Astrophys. 
J. 517 (1999) 565; W. Michael Wood-Vasey et al., Astrophys. J. 666 (2007) 
694; Pierre Astier et al., Astron. Astrophys. 447 (2006) 31. 

[2] E. Komatsu et al., larXiv:0803.0547l [astro-ph]. 

[3] M. Tegmark et al., Astrophys. J. 606 (2004) 702; M. Tegmark et al., Phys. 
Rev. D 69 (2004) 103501. 

[4] T. Padmanabhan, Phys. Rept. 380 (2003) 235; E.J. Copeland, M. Sami, S. 
Tsujikawa, Int. J. Mod. Phys. D 15 (2006) 1753; J. Frieman, M. Turner, D. 
Huterer. la7Xiv:0803.0982l [astro-ph]: T. Padmanabhan, arXiv:0807.2356 [gr- 
qc]. 

[5] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1. 

[6] N. Straumann, in Duplantier, B. (ed.) et al.: Vacuum energy, renormalization, 



7-51, arXiv:astro-ph/0203330 



[7] S. Nobbenhuis, arXiv:gr-qc/0609011 



[8] H. Stefancic. larXiv:0807.36"92l [gr-ac], to appear in Phys. Lett. B. 

[9] S. Nojiri, S.D. Odintsov, Phys. Rev. D 72 (2005) 023003. 

[10] S. Capozziello, V.F. Cardone, E. Elizalde, S. Nojiri, S.D. Odintsov, Phys. Rev. 
D73 (2006) 043512. 

[11] I. Brevik, E. Elizalde, O. Gorbunova, A.V. Timoshkin, Eur. Phys. J. C52 
(2007) 223. 

[12] I. Brevik, O.G. Gorbunova, A.V. Timoshkin, Eur. Phys. J. C51 (2007) 179. 
[13] K. Bamba, S. Nojiri, S.D. Odintsov, JCAP 0810 (2008) 045. 



11 



[14] D.F. Mota, C. van de Bruck, Astron. Astrophys. 421 (2004) 71. 

[15] S. Nojiri, S.D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4 (2007) 115, 
|hep-th/0601213| 

[16] T.P. Sotiriou, V. Faraoni, i?Xm0805.1726 [gr-qc]. 

[17] S. Nojiri, S.D. Odintsov, Phys. Rev. D 68 (2003) 123512. 

[18] G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, L. Sebastiani, S. Zerbini, 
Phys. Rev. D77 (2008) 046009. 

[19] S. Weinberg, Astrophys. J. 168 (1971) 175. 

[20] W. Zimdahl, Phys. Rev. D 53 (1996) 5483. 

[21] 0. Gr0n, Astrophys. Space Sci. 173 (1990) 191. 

[22] R. Colistete, Jr., J.C. Fabris, J. Tossa, W. Zimdahl, Phys. Rev. D76 (2007) 
103516. 

[23] A. Avelino, U. Nucamendi, F.S. Guzman, AIP Conf. Proc. 1026 (2008) 300. 

[24] A. Avelino, U. Nucamendi. larXiv:0810^0303l [gr-qc], 

[25] I. Brevik, O. Gorbunova, Gen. Rel. Grav. 37 (2005) 2039. 

[26] I. Brevik, S.D. Odintsov, Phys. Rev. D65 (2002) 067302. 

[27] I. Brevik, Int. J. Mod. Phys. A18 (2003) 2145. 

[28] I. Brevik, O. Gorbunova, Y.A. Shaido, Int. J. Mod. Phys. D14 (2005) 1899. 

[29] I. Brevik, J.-M. Borven, S. Ng, Gen. Rel. Grav. 38 (2006) 907. 

[30] I. Brevik, Gen. Rel. Grav. 38 (2006) 1317. 

[31] I. Brevik, Int. J. Mod. Phys. D15 (2006) 767. 

[32] I. Brevik, Eur. Phys. J. C56 (2008) 579. 

[33] I. Brevik, larXiv:0808.3644l [gr-qc] 



12 



